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ABSTRACT 


The  advent  of  tilt  rotor  technology  asks  for  rotors  that 
have  different  twist  and  RPM  requirements  in  hover  and  in 
forward  flight  to  optimize  for  operational  conditions.  In 
order  to  get  an  assessment  of  the  capabilities  to  fulfill 
these  requirements  this  report  presents  a  mapping  of  twist 
angle  variation  as  a  function  of  RPM  and  laminate  orientation. 
The  basic  laminate  for  the  six  models  as  well  the  D-shape  spar 
that  represents  the  structurally  active  part  of  the  blade  is 
assumed  to  be  constant  ( 0°/90o/0o/90o/90°/0o ) .  This  six  layer 
cross  ply  laminate  is  chosen  as  it  provides  the  necessary 
extension-twist  coupling  without  a  hygrothermally  induced 
twist  that  is  highly  undesirable.  The  couplings  and  trends  in 
the  models  are  visualized  in  carpet  plots,  one  for  each  model, 
in  an  attempt  to  establish  a  method  to  answer  the  basic 
question  of  the  magnitude  of  twist  angle  available  due  to  a 
particular  geometry,  material  and  load  system. 
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I NTRODUCT I ON 


A.  .  OVERVIEW 

In  the  design  process  of  composite  material  structures  it 
is  often  noticed  that  the  best  characteristics  of  the  material 
are  not  used.  The  normal  procedure  is  to  tailor  the  material 
to  approach  the  properties  of  isotropic  materials,  avoiding 
the  "undesirable"  couplings  between  extension  and  twist  for 
example . 

That  traditional  approach  worked  for  helicopter  design. 
However,  in  the  tilt  rotor,  the  idea  of  changing  the  twist 
angle  of  the  blade  to  optimize  for  either  hover  or  forward 
flight  came  into  the  scenario. 

This  study  will  use  a  blade  spar  to  study  the  sensitivity 
of  the  extension- twist  coupling  due  to  changes  in  RPM,  layer 
angles  definition,  changes  in  laminate  orientation  and  also 
the  effects  of  each  of  these  configurations  in  the  first  four 
natural  vibration  modes  of  the  blade. 

In  the  analysis,  a  program  named  CASA/GIFTS  (Computer 
Aided  Structural  Ana lysis /Graphical  Interactive  Finite  Element 
Total  System)  is  adopted.  The  program  makes  it  easier  to 
visualize  the  effects  of  the  coupling  in  the  structure  and  is 
seen  to  be  an  efficient  design  tool.  The  results  obtained 
through  the  finite  element  procedure  for  each  laminate 
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configuration,  in  a  total  of  six  models  will  be  presented  in 
3-D  graphs  that  show  the  trends  for  each  model  for  the  layer 
angle  (top  layer  as  a  reference),  the  RPM  and  the  twi3t  angle. 

These  graphs  can  then  be  used  as  a  preliminary  orientation 
by  the  designer  to  achieve  the  optimal  or  near  optimal 
combination  of  laminate  configuration  layer  angle  within  the 
laminate  that  fit  the  specific  requirements  for  the  particular 
blade  in  focus. 

B  .  literature:  review 

This  section  reviews  briefly  some  pertinent  publications 
listed  as  References  in  the  thesis  report.  Other  related 
literature  found  during  the  scope  of  this  research,  is  listed 
under  Bibliography.  The  topics  reviewed  are  classified  under 
the  following  categories: 

.  Finite  Element  Theory. 

.  Composite  Material  Theory. 

.  Vibrations. 

.  Helicopter  Theory. 

.  Design  Applications  (Modeling). 

.  Auxiliary  Software. 

In  the  Finite  Element  Theory,  relevant  to  present 
research,  some  basic  references  are:  Batozfl980]  that 
describes  the  finite  element  procedures  in  engineering 
analysis  using  Discrete  Kirchoff  Theory  ( DKT ) ,  Cookf 19811  on 
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basic  concepts  and  applications  of  finite  element  analysis, 
Craig[1981]  describing  computer  methods  in  structural  dynamics 
and  f CASA/GIFTS,  1987]  that  is  the  reference  manual  for  the 
finite  element  code  used  for  the  present  (static  and  dynamic) 
work.  Its  graphical  capabilities  are  very  convenient  to  the 
analyst  in  a  manner  that  simplifies  the  visualization  of  the 
effects  of  loads  and  couplings  in  the  structure,  making 
difficult  concepts  ruore  accessible  to  understanding. 

Tsai  and  Pagano[1968]  establish  a  notation  in  which  the 
composite  lamina  properties  are  invariant  with  respect  to  the 
axis  of  rotation.  This  approach  provides  a  very  useful  way  to 
compare  various  material  systems;  the  laminate  theory  is  well 
documented,  see  for  example  Vinsonfl987]  and  Jones[1975]  both 
offer  a  good  source  on  the  mechanics  and  behavior  of  composite 
materials . 

Yntema[1955]  developed  a  very  useful  tool  to  estimate 
bending  frequencies  of  rotating  beams  in  which  the  stiffness 
effect  of  RPM  is  shown  for  several  blades  (beam) 
configurations.  Schilhansl T 1958  ]  uses  another  approach  to 
estimate  bending  frequencies  of  rotating  cantilever  beams, 
introducing  the  effect  of  the  angle  made  by  the  minor  axis  of 
inertia  with  the  direction  of  the  tangential  circular 
velocity.  Wood[  19651  shows  a  parametric  investigation  of 
aerodynamic  and  aeroelastic  characteristics  of  articulated  and 
hingeless  rotor  systems.  Pritchardf 19881  presents  an  optimal 
placement  of  tuning  masses  to  reduce  vibration  levels  in 
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helicopter  blades.  Kottapalli [ 1983 1  approaches  the  vibration 
reduction  problem  by  modifying  the  blade  torsional  response. 
This  is  an  innovative  approach  in  that  he  modifies  the 
aerodynamic  response  instead  of  adding  mass,  which  is  the 
usual  conventional  way  to  solve  the  vibration  reduction 
problem. 

The  helicopter  theory  in  general  is  found  comprehensively 
in  Johnsonf 1980 3  and  is  complemented  in  this  research  by 
Prouty[1986]  in  the  understanding  design  trends. 

Hoskin[1986]  is  a  general  informative  book  in  the 
aeronautical  utilization  of  composite  materials,  and  in 
Lake[1988]  there  is  a  preliminary  investigation  of  finite 
element  modeling  of  composite  material  rotor  blades. 
Nixon[1987]  addresses  the  extension-twist  coupling  of 
composite  circular  tubes  applied  to  tilt  rotor  design.  In 
McVeigh  F 1983]  the  aerodynamic  design  of  XV-15  advanced 
composite  tilt  rotor  blade  requirements  are  reviewed.  Finally 
Hodgesfl987]  presents  a  comparison  of  composite  material  rotor 
blade  models,  using  two  different  methods  of  analysis:  this 
model  was  chosen  as  a  prime  reference  to  build  the  basic  model 
used  throughout  this  work. 

Several  software  tools  made  it  possible  to  manipulate  the 
immense  volume  of  data  transforming,  it  into  a  suitable  and 
compact  presentation  form,  enabling  the  conclusions  to  be 
obtained.  Kellyf 19883  was  the  basic  reference  to  use  Word 
Perfect  5.0  in  the  word  processing  and  graphical  editing  jobs. 
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Holt [19881  is  a  reference  that  helps  to  manage  the  data  '  .vto 
Lotus  1-2-3  files  and  to  get  the  input  graph  files  to  Word 
Perfect  and  to  SURFERf 1 987 ] .  The  SURFER  program  was  used  to 
obtain  the  3-D  carpet  plots  of  extension-twist  coupling 
analysis . 

C.  THESIS  OUTLINE 

This  thesis  report  is  divided  into  three  areas,  often 
found  during  a  design  process. 

It  begins  by  modeling  a  rotor  blade  with  laminate 
composite  material.  A  numerical  solution  is  obtained, 
exploiting  the  extension-twist  coupling  inherent  when 
asymmetries  are  present  in  the  laminate  construction. 
Discrepancies  encountered  in  the  results  compared  with  other 
rotor  models  are  reported.  These  differences  raised  questions 
fostering  the  research  towards  the  possibilities  of  different 
rotor  constructions,  starting  with  the  same  basic  laminate. 

The  different  rotor  configurations  obtained  are  then 
subject  to  both  static  and  dynamic  analysis.  The  static 
analysis  reveals  the  different  responses  of  these  models, 
while  the  dynamic  behavior  remains  almost  invariant,  yielding 
only  small  variations  in  the  first  four  fundamental  vibration 
modes . 

Analysis  of  the  results  is  done,  and  it  is  believed  that 
a  new  insight  into  finite  element  analysis  of  composite 
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laminates  is  provided.  This  new  insight  is  better  noticed  when 
there  are  asymmetries  present  within  the  laminate 
construction. 
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1 1 


THEORETICAL  FORMULATION 


A.  BASIC  CONSIDERATIONS 

To  understand  how  composite  material  works  and  its 
particular  behavior  that  is  useful  in  "designing"  the 
material,  a  macroscopic  approach  is  re  sorted  to. 

Composite  means  that  two  or  more  materials  are  combined 
and  usually  preserve  their  best  qualities,  and,  often  exhibit 
properties  that  neither  one  possesses  by  itself. 

Some  of  the  properties  that  can  be  improved  by  making  a 
composite  material  are: 

.  strength 
.  stiffness 

.  corrosion  resistance 
.  wear  resistance 
.  life  in  fatigue 
.  thermal  insulation 
.  thermal  conductivity 
.  acoustical  insulation 

It  may  be  noted  that  not  all  these  properties  are  improved 
at  the  same  time,  and  usually  it  may  not  be  required. 
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An  Orthotropic1  material  has  properties  which  are 
different  in  three  mutually  perpendicular  directions  at  a 
point,  in  other  words  has  three  different  planes  of  symmetry. 

A  lamina  is  the  basic  form  of  any  laminated  fiber- 
reinforced  composite,  it  consists  of  a  plane  mounting  of 
fibers  imbedded  in  a  matrix. 

In  the  macro-mechanical  behavior  of  the  lamina,  only 
averaged  apparent  mechanical  properties,  or  better,  the 
stress-str? in  relations  for  an  orthotropic  material  under 
plane  stress  conditions,  transformed  to  directions  not  aligned 
with  the  principal  directions  of  the  lamina,  are  considered. 

This  describes  the  behavior  of  the  orthotropic  material 
whose  laminae  have  arbitrary  directions,  other  than  the 
natural  geometric  directions  of  the  structure  to  be  analyzed. 

13  .  STRESS  — STRAIN  RELATIONS 


The  stress-strain  relations  or  generalized  Hooke’s  law  for 
anisotropic  material  will  be  reviewed  briefly,  and  the  basic 
restriction  of  linear  elastic  behavior  assumption  applies. 

The  following  derivations  can  be  found  in  more  detail  in 
Vinson[ 1987 ]  . 

The  elasticity  tensor  has  nine  independent  elasticity 
constants,  because  of  symmetry  (Cij=Cji)  and  the  three  mutually 


Orthogonally  anisotropic. 
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orthogonal  planes  of  symmetry  (C16=C26=C36=C45=0 ) ; 
that  is, 


Cu 

c12 

c13 

0 

0 

0 

C21 

c22 

C23 

0 

0 

0 

C31 

C32 

C33 

0 

0 

0 

0 

0 

0 

C<4 

0 

0 

0 

0 

0 

0 

C55 

0 

0 

0 

0 

0 

0 

C, 

and 


(2) 


*i  =  x  6j  (  i  =  l ,  2 , 3 , 4 , 5 , 6  ;  j  =  l ,  2 , 3 , 4 , 5 , 6  )2 


In  these  relations  the  hygrothermal  and  thermal  effects 
are  not  included  for  the  sake  of  simplicity  of  this  report. 

Through  a  force  equilibrium  study,  analogous  to  the  Mohr’s 
circle  analysis  [Jones,  1975]  the  stresses  and  strains  in  the 


2The  material  directions  1  to  6  are  equivalent  to 
cartesian  directions:  x,y,z,yz,xz,xy  respectively  in  material 
axis . 
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material  directions  can  be  related  to  arbitrary3  directions 
using  the  following  relations: 


•  “ 

” 

°x 

CTl 

°y 

°z 

=  [T]-X 

o2 

°3 

°yz 

°xz 

°5 

°x  y 

°6 
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“ 

ey 

€2 

6z 

=  [T]'1 

e3 

* yz 

*4 

^xz 

e5 

exy 

66 

The  transformation  matrix  is  constructed  by  defining 
m=cos  0,  n=sin  0;  where  0  is  the  ply  orientation  angle, 
defined  with  respect  to  material  and  structural  axes,  assuming 
positive  sign  in  counterclockwise  direction  (Figure  1). 


3  More  convenient,  xyz  axis  related  to  the  structure 
geometry  for  example. 
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Figure  1  Positive  rotation  of  principal  material  axes  from 
arbitrary  xy  axes  [Jones,  1975] . 


The  inverse  of  transformation  matrix  [T]  is  expressed  in 
terms  of  direction  cosines  m  and  n  as: 
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For  a  lamina  of  composite  material  the  modern  notation 
uses  for  the  quantities  in  the  material  stiffness  matrix 
and  are  defined  as: 


(6) 


11s 

=E11 

(1~v23v32 

)/6, 

®22=E22*1' 

v31v13)/6 

33  = 

=E33 

(1-v12v21 

)/8, 

®44=  G23’ 

Q55 

=  G13 

12: 

=E11 

(v21+v31v 

23  )/6 

*  E22(v12+v 

32V13 

)/6 

13: 

=EU 

(v31+v21v 

32  >/6 

“  E22(vX3+v 

12V23 

)/8 

23  = 

=E22 

<v32+v12v 

31)/6 

=  E33(v23+v 

21V13 

)/6 

5  -  1  -  v12v21-  v23v32-  v31v13-  2v21v32v13 


If  the  lamina  has  same  properties  in  both  2  and  3 
directions,  then  viz=vi3,  v31=v32=0,  G12=G13,  E22=E33  and  a 
simplification  in  equations  (6)  lead  to  equations  (7). 


Q11=E11/(1  "  v12v21)’  Q22=E22/(1  "  V12V21) 


(7)  Q12=  v2 i ei l ^ ^ 1  "  v12v21^=  v12E22  ~  V12V21^ 


Q66=  G12 
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To  compare  various  material  systems  in  the  design  of  a 
composite  structure,  Tsai  and  Pagano[1969]  arrived  at  some 
invariants  with  respect  to  the  axis  of  rotation: 


V  1  <3®n+3«22+2®12+4«66> 


V 

1 

2 

<«»- 

Q22> 

(8) 

V 

1 

8 

(®n+ 

Q22“ 

2Q12- 

4Q66 

V 

1 

8 

(«11- 

Q22  + 

6Q12- 

4Q66 

V 

1 

8 

(®11+ 

Q22“ 

2*12  + 

4Q66 

These  invariants  are  representative  of  the  lamina 
properties.  The  [Q]  matrix  may  be  transformed  to  x-y  axes  from 
1-2  axes  by  the  usual  coordinate  transformation. 


(9) 

[Q] 

=  IT]"1  [Q][T] 

In  terms  of 

the 

invariants 

are 

written  as: 

®11  = 

V 

UgCosf  20 ) 

^22  = 

V 

UgCos ( 20 ) 
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After  the  definitions  in  (10),  the  stress-strain  relations 
for  the  kth  lamina  of  a  N  ply  laminate  (Figure  2)  are  given  by: 
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Where  are  defined  earlier. 


Figure  2  Laminate  construction,  [ 90°/+45o/0o/-45o/90° ] 
characterized  in  this  case  to  show  the  nomenclature  for 
the  lay-up. 
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LAM INATE  ANALYSIS 


When  two  or  more  laminae  are  bonded  together  to  obtain  an 
integral  structural  component,  the  laminae  principal  material 
directions  are  established  to  produce  a  component  that  can 
resist  loads  in  several  directions  (Figure  2). 

The  following  procedures  enable  to  obtain  the  stiffness 
matrix  of  such  a  composite  material  from  the  lamina 
properties,  using  arbitrary  combinations  of  materials  and 
layer  angles. 

For  the  kth  lamina  of  a  N  layer  laminate,  equation  (11) 
represents  the  stress-strain  relations. 

In  order  to  obtain  the  strain-displacement  relation.  '•.he 
displacements  are  defined  as  follows: 


(12) 


u(x,y,z ) 


v(x,y,z) 

w(x,y,z) 


u0(x,y)  +  za(x,y) 
v0(x,y)  +  z(3(x,y) 
w(x,y) 


where  u,  v  and  w  are  displacements  in  x,  y  and  z  directions 
respectively 
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and4 


6w 


6w 


(13)  a  =  -  —  and  0  =  -  — 


6x 


6y 


are  the  negative  of  the  slope,  related  to  the  rotations. 


The  mid-surface  strains,  then  are  given  by 
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while  the  curvatures  are  given  by 
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4  In  the  following  expressions  6  means  partial 
derivative , 
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Defining  h  as  laminate  thickness,  hk  (Figure  3)  as  the 
vectorial  distance  from  the  panel  mid-plane,  N,  M,  Q  as  stress 
resultants  (Figure  4),  stress  couples  and  shear  resultants 
respectively.  Relations  can  be  established  relating  these 
resultants  to  the  strains  and  curvatures.  This  matrix  is 
henceforth  referred  as  the  ABD  matrix. 

This  nomenclature  is  adopted  for  clarity  and  simplicity 
as  may  be  noted  from  the  expressions  that  follow. 

The  stress  resultants  are  expressed  by  integrating 
stresses  across  the  thicicness  of  the  plate  and  is  given  by: 
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and  the  moment  resultants  take  the  following  form 
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Figure  3  Nomenclature  for  the  Stacking  Sequence 
[Vinson,  1987  ] . 


Figure  4  Positive  directions  for  Stress  Resultants  and 
Stress  Couples  for  a  Plate  [Vinson, 1987] . 
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Relating  resultants  with  the  inplane  strains  and 
curvatures  through  expressions  111)  through  (14),  matrix  form 
is  given  by: 


IN]  =  [A].[e0]  +  [B].[K] 

(18) 

[M]  =  [B] . [e0]  +  [D] . [K] 


In  expression  (18),  the  matrices  A,  B  and  D  are  called 
extensional  stiffness  matrix,  bending-stretching  coupling 
matrix  and  flexural  stiffness  matrix  respectively,  where: 


Aij 

=  (  ®ij>k ' 

hk  '  hk-l’ 

[i, j  =  1 , 2 , 6 ] 

(19) 
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Dij 

■  1  *  2k=i(  ®ij 

>k  1  hk  -  hk-l’ 

[i, J  =  1 , 2 , 6 ] 

Assuming  that  transverse  shear  stresses  are  distributed 
parabolically  (even  though  there  exists  discontinuities  at 
the  laminae  interfaces),  the  shear  resultants  may  be  written 
as : 
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where 


(21) 
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®ij>k 


i,j  =  4,5  only 

The  stresses  and  moments  resultants  can  be  written  in  a 
compact  form  that  is  very  convenient  to  see  the  overall 
relations  within  a  composite  material  laminate,  specifically 
the  ABD  matrix.  The  ABD  matrix  is  obtained  as  a  combination 
of  the  expressions  in  (18).  Expressions  (18)  and  (22)  are 
equivalent  and  in  expression  (22)  the  ABD  matrix  is  defined 
explicitly. 
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D -  EFFECTS  OF  FLY  ANGLE  © 


To  get  an  insight  of  the  overall  influence  of  the  ply 
angle  8  on  the  laminate  stiffness,  it  is  useful  to  observe 
the  effects  of  0  for  a  given  lamina.  As  each  element  of  the 
ABD  matrix,  Ai{ ,  Bi;j,  and  Dy ,  is  a  function  of  the  material 
stiffness  coefficient  Qii?  the  layer  thickness  tk  and  mid-plane 
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vectorial  distance  hk,  the  Q^’s  represent  the  stiffness 
variations  due  to  ply  angle  changes  in  a  given  lamina. 

The  property  of  the  composite  material  that  reveals  its 
advantages  over  metallic  or  plastic  is  that  it  can  be  tailored 
to  fit  a  given  structural  shape  and  a  specific  set  of  loads 
and  operating  conditions. 

The  terms  of  the  ABD  matrix  are  functions  of  the  material 
properties  and  the  ply  angle  0.  If  the  properties  are  kept 
constant  and  within  typical  values  of  most  fiber  resin  types 
of  composite,  it  is  possible  to  plot  the  values  of  versus 
0  and  get  a  set  of  seven  very  useful  graphs.  These  will  permit 
an  initial  appreciation  of  such  effects,  and  how  they  can  be 
combined  to  meet  a  specific  structural  requirement. 

Figures  5  to  11  show  how  the  elements  of  the  "material 
stiffness"  matrix  ABD  can  be  tailored  by  changing  the  ply 
angle,  taking  full  advantage  of  this  characteristic  behavior 
of  laminated  composite  material. 
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In  figure  11  all  Q’s  are  superimposed,  on  the  same  plot, 
to  give  a  better  perspective  of  the  trends  in  doing  this  kind 
of  layer  design  to  fit  specific  structural  constraints.  When 
A16  and  A26  are  non-zero  there  is  stretching-shearing  coupling, 
while  non-zero  B16  means  twisting-stretching  coupling,  and  non¬ 
zero  B26  indicates  bending-shearing  coupling;  bending-twisting 
coupling  is  originated  in  non-zero  values  of  D16  and  D26. 

In  normal  design  procedures,  these  terms  are  made  zeroes 
by  selecting  appropriate  stacking  sequence;  cross  ply  (0°/90° 
plies  combination)  or  angle  plv  laminates  (-G/+0  plies 
combination).  Other  important  factor  to  be  considered  is  the 
symmetry  with  respect  to  the  mid-surface  plane. 
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The  stress  strain  relations  or  constitutive  relations  (ABD 
matrix),  combined  with  the  proper  equations  of  equilibrium  and 
the  strain-displacement  relations  form  the  basis  for  analyzing 
beams,  plates  or  shells. 

These  basic  relations  lead  to  the  motivation  and  reasoning 
behind  the  model  tailoring  to  fit  a  design  requirement. 


Ill 


ROTOR.  MODELING 


A.  TYPI CAL  CURRENT  ROTOR 

MODELING 


Current  rotor  design  is  a  compromise  between  the  different 
requirements  in  hover  and  forward  flight,  and  the  analytical 
models  reflect  that  tendency  or  limitation.  For  instance,  the 
optimum  twist  for  hover  is  not  the  same  in  forward  flight  for 
tilt  rotors  as  well  as  for  helicopters.  Indeed,  very  often, 
the  designers  avoid  extension- twist  coupling  with  laminated 
composite  materials  because  of  certain  undesirable  effects. 
The  result  is  that  one  of  the  two  flight  modes  will  have  less 
than  optimal  propulsive  efficiency,  depending  on  design  phase 
decisions  driven  by  requirements  and/or  aircraft  mission 
profiles . 

With  the  tilt  rotor  concept  of  mixing  fixed  wing  and 
rotorcraft  technology,  these  design  trends  were  brought  to 
surface  once  more. 

B.  PROPOSED  SOLUTION 

Nixon[1987]  proposed  a  solution  for  this  sort  of  dilemma. 
The  model  used  an  extension-twist  coupled  rotor  blade  with 
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100%  RPM5  in  hover  and  80%  RPM  in  forward  flight.  These  are 
reasonable  ratios  between  hover  and  forward  flight  for  a  tilt 
rotor  type  aircraft.  The  difference  in  RPM  provides  a  change 
in  centrifugal  force,  inducing  a  change  in  the  twist  angle 
distribution,  which  can  be  used  to  optimize  the  rotor  for  both 
regimes . 

The  use  of  composite  materials  makes  it  possible  to  design 
a  rotor  blade  that  takes  advantage  of  the  changing  loads, 
rendering  the  required  twist  distribution  optimum  in  each 
flight  mode.  At  this  point  one  question  arises: 

What  is  the  magnitude  of  twist  available  from  a  particular 

geometry  and  material  system? 

The  answer  to  that  question  lies  in  a  method  of 
prediction,  with  a  model  that  behaves  as  close  to  the  actual 
blade  as  possible  yet  inexpensive  enough  to  permit  a  wide 
spectrum  parametric  study,  providing  a  source  to  evaluate 
several  geometry  and  material  systems.  Certainly  such  modeling 
should  be  verified  by  experimental  data. 

This  research  is  aimed  to  address  the  analytical  part  of 
the  above  question  and  provide  design  tools  to  rotor  blade 
systems . 


5  Reference  value  for  100%  RPM  is  217  rad/s,  in  that 
scale  of  rotor  blade. 
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The  torsional  moments  that  affect  a  rotating  blade  can 
modify  substantially  the  twist  distribution. 

In  this  numeric  model,  the  pitching  moments  generated  by 
aerodynamic  loading  are  not  taken  into  account,  the  emphasis 
being  on  the  inertia  and  extension-twist  coupling  as  sources 
of  torsional  moments. 

An  assumed  lift  distribution  acts  on  the  model,  but 
compared  with  the  inertia  forces,  its  effects  are  negligible 
on  the  rotating  blade.  The  inertia  forces  acting  on  the 
rotating  blade  produce  extension  which  causes  twist  by 
structural  coupling.  However,  the  same  inertia  forces  produce 
a  centrifugal  flattening  (Prouty,  1986]  effect  called  tennis 
racket  effect  (because  of  the  tendency  of  a  tennis  racket  to 
align  its  plane  with  the  plane  of  rotation  as  it  is  swung  in 
an  arc ) . 

The  forces  acting  in  a  helicopter  blade  can  be  seen  in 
Figure  12,  and  it  may  be  noted  that  in  this  case  the  effect 
is  reduced  with  balancing  weights  (Chinese  weights). 

O.  GIFTS  CAPABILITIES 

Gifts  capabilities  include  generation  of  any  type  of 
structural  model,  with  the  following  elements: 

.  one-dimension  rods  or  beams 

.  two-dimension  plates 
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.  three-dimension  shells,  solid  and  complex  stiffened 

shells . 

The  elements  can  be  selected  from  a  library  of  options. 

The  materials  can  be  created  by  the  user  or  selected  from 
a  library  of  standard  definitions.  The  program  allows  the  user 
to  define  isotropic  materials  and  laminated  composite 
materials,  that  can  have  up  to  a  maximum  of  hundred  layers. 

Loads  can  be  concentrated  or  distributed,  static  as  well 
as  function  of  time. 

Static  and  dynamic  analysis  are  performed  on  the  models 
generated.  The  dynamic  analysis  provides  free  vibration 
response  and  forced  response  (transient  computations),  and 
model  superposition. 

In  using  GIFTS  as  the  finite  element  program  to  model  the 
blade,  the  effects  of  the  centrifugal  force  are  taken  into 
account  with  the  forces  acting  on  the  lumped  masses  as  shown 
in  Figure  12,  radially  from  the  center  of  rotation.  In  using 
other  finite  element  programs  this  feature  must  be  verified 
for  correct  usage  and  application.  The  tennis  racket  effect 
can  affect  the  twist  angle  to  a  large  extent.  If  the 
centrifugal  forces  act  in  parallel  lines  instead  of  passing 
through  the  center  of  rotation,  there  will  be  no  centrifugal 
flattening  effect  and  results  in  larger  twist  angles. 

The  basic  rotor  model  selected  to  apply  the  proposed 
method  is  the  one  analyzed  by  Hodgesf 1987 1 ,  this  model  was 
analyzed  using  two  different  methods,  namely  by  finite  element 
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( MSC/NASTRAN )  and  a 
coupled-beam  model 
[Rehfield,  1985].  The 
rotor  geometry 
presented  by 
Hodges[1987]  was 
modeled  in  GIFTS  using 
triangular  TB3  element 
( "  flat  triangular 
plate  or  shell 
element,  including 
membrane,  out  of  plane 
and  fictitious  in- 
plane  bending 
stiffness . " [Batoz , 
1980]).  The  results 
indicated  a  wide 
disparity  of  results. 
However,  the  present 
model  appears  to  be 
within  the  assumption  of  linear  and  small  displacements 
theory . 

The  models  in  Figure  13  have  a  notation  (M120  /  Ml-20) 
that  is  described  in  the  next  section.  A  twist  angle 
distribution  discrepancy  could  be  explained  by  the  tennis 
racket  effect  probably  not  present  in  the  model  described  in 


Figure  12  Twisting  Moments  Due  to 
Centrifugal  Forces  fProuty, 1986 ] . 
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Hodgesf 1987 ] ,  but  the  tip  vertical  displacement  due  to  lift 
and  weight  with  no  rotation  could  not  be  explained  by  this 
argument.  It  may  be  observed  that  the  displacements  quoted  in 
the  referenced  work  are  about  25%  of  the  length,  which  appears 
to  violate  the  linear  theory  that  was  used. 
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BLADE  .DEFLECTION  DUE  TO  LIFT  AND  WEIGHT 


COMPARISON  WITH  NASA  TM  8802* 


Figure  13  Comparison  between  GIFTS  model  and  NASA  TM  89024 
resul ts. 


This  discrepancy  may  be  due  to  a  possibility  that  the 
models  compared  are  not  exactly  the  same,  although  every 
effort  was  made  to  take  into  account  all  the  pertinent 
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geometric  and  material  data  reported.  However  this  discrepancy 
was  the  driving  factor  that  raised  some  questions,  such  as: 
What  are  the  possibilities  to  model  different  structural 
configurations  starting  with  the  same  number  of  layers, 
same  lay-up  design,  same  thickness  and  material 
properties?  (perhaps  by  inadvertent  oversight  by  the 
designer,  analyst  using  canned  programs,  or  manufacturing 
engineer.  ) 

The  answer  to  this  was  based  on  the  following  assumptions: 

.  The  structure  is  a  D  shape  spar. 

.  It  is  built  using  three  shells,  top,  bottom  and  rear 
shell  surfaces  . 

.  The  top  and  bottom  meet  at  the  leading  edge  and  the  rear 
shell  closes  the  box. 

.  Once  established,  the  internal  lay-up  in  the  stacking 
sequence  remains  constant. 

With  those  assumptions  a  laminate  to  constitute  one  shell 
surface  may  be  defined,  using  only  the  Laminate  Principal 
Orientation  ( LPO )  angle.  For  example,  in  the  following 
laminate  L  +  20°/-70o/  +  20o/-70o/-70o/+20°]  ,  the  top  layer  angle 
is  denoted  +20°  and  is  a  six  layer  asymmetric  laminate. 

Rotating  the  whole  laminate  40°  clockwise  gives 
[ -20°/  +  70o/-20o/  +  70o/  +  70o/-20°]  ,  whose  Laminate  Principal 
Orientation  angle  is  denoted  as  -20°. 

This  laminate  has  intentionally  designed  asymmetry  to  get 
the  ex tension- twi st  coupling.  Associated  with  this  laminate 
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is  a  normal  vector  whose  orientation  is  seen  in  Figure  14.  The 
D-spar  then  is  built  using  three  such  laminates,  with  each 
laminate  mounted  with  the  normal  vector  facing  inward  or 
outward.  This  method  of  construction  results  in  six  different 
configurations  of  the  rotor.  These  configurations  and  other 
variations  within  each  one  will  be  developed  in  the  next 
section . 

D.  VARIOUS  ROTOR  DESIGNS  OR 
CONFIGURATIONS 


The  six  different  configurations  are  designated  Ml  to  M6 
and  the  LPO6,  varying  from  -90°  to  +90°  in  steps  of  20°, 
provides  the  suffix  to  complete  the  model  designation  Ml-907, 
M340,  M6-60  etc.  In  a  total  of  11  LPO  in  each  model. 

Different  configurations  may  be  visualized  with  the  help 
of  the  scheme  shown  in  Figure  14. 

Model  I  and  IV  can  be  built  with  only  one  laminate  wrapped 
around,  one  inverse  of  the  other. 

In  the  program  GIFTS,  the  outward  vector  for  each  surface 
is  defined  during  the  "GRID"  definition,  e.g.  in  a  four  sided 
grid,  each  shell  is  defined  within  a  grid.  The  local  X  axis 


6  Laminate  Principal  Orientation  Angle,  referred  to  the 
top  layer  of  the  laminate  used  to  construct  that  particular 
model . 

7  Ml-90  meaning  model  1  with  LPO  =  -90°;  M340  meaning 
model  3  with  LPO  =  40°. 
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is  oriented  parallel  to  the  first  line  of  the  grid,  the  local 
Z  axis  is  the  cross  product  of  the  first  and  second  "line"  in 
the  grid  definition,  while  the  local  Y  axis  is  the  cross 
product  of  the  Z  and  X  axes,  as  is  usual  in  cartesian  systems. 

The  outward  vector  is  oriented  in  a  direction  defined  by 
the  local  Z  axis.  Layer  angles  are  positive  in  the  counter¬ 
clockwise  direction  in  the  XY  plane  for  shell  elements,  as  is 
the  convention  for  laminated  composite  materials. 

END  VIEW  OF  D  SPAR  (  3  SHELLS) 

r  - ^  - 

— - - 

M0DELX  1  MODEL  2  MODEL  3 

-E3^ 

MODEL  4  MODEL  5  MODEL  6 
NORMAL  VECTOR  ASSOCIATED  TO 
LAMINATE  STACKING  SEOUENCE 

Figure  14  Normal  to  surface  vector  associated  with  the 
layer  sequence  of  the  basic  laminate  and  its  position  in 
z'sJ.&’tion  to  end  view  of  the  bi sde • 

It  may  be  observed  that,  when  an  asymmetry  is  present  in 
the  basic  laminate  (material  definition),  the  structural 
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response  is  highly  dependent  on  how  the  grid  and  consequently 
the  lay-up  is  defined.  In  the  case  of  a  symmetric  laminate  or 
isotropic  material,  this  fact  is  never  brought  to  surface 
because  of  its  irrelevancy.  However,  dealing  with  asymmetric 
laminated  composite  materials,  the  analyst  must  always  have 
this  in  mind.  It  may  be  noted  that,  in  most  cases  the  laminate 
is  chosen  to  be  symmetric. 

Table  1  gives  in  a  compact  form  the  description  of  various 
laminate  constructions  described.  Included  in  that  table  are 
the  line  names  actually  employed  in  the  input  file  M$**.SRC 
[Annex  1]  ($  =  1,2, 3, 4, 5, 6  models;  **  =  top  layer  angle  as 
defined  previously,  in  the  range  -90  to  +90),  to  develop  the 
models  generated  in  GIFTS.  By  modifying  two  areas  in  this 
input  file,  six  models  are  generated  that  are  considered  under 
the  present  study,  this  also  provides  variation  of  the  layer 
angle  within  each  model.  The  actual  modifications  employed  in 
the  input  file  are  given  in  LAYUP. AUX  [Annex  2].  For  each 
model  (one  of  the  six  combinations  obtainable  by  the  change 
in  the  laminate  associated  vector,  as  given  in  Figure  14), 
the  whole  laminate  was  rotated  in  increments  of  10°, 
referenced  to  the  top  layer  angle  (LPO).  For  each  of  these 
models ,  the  RPM,  and  consequently  the  centrifugal  force,  was 
varied  in  the  range  (130  rad/s  to  300  rad/s). 

The  model  geometric  characteristics  are  the  same  as  in 
Hodgest 1987  ]  ,  a  graphite  epoxy  composite  rotor  blade 
constructed  out  of  Hercules  IM6  fiber  with  Ciba-Geigy  R6376 
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resin.  The  cured  ply  thickness  is  assumed  to  be  0.0055  in.  The 
D-spar  has  a  35.23  in.  radius,  and  begins  at  station  5.32  in., 
which  has  all  degrees  of  freedom  suppressed  to  emulate  a  rigid 
rotor  connection  to  the  hub.  The  rotor  has  a  constant  cross- 
section,  defined  by  13  nodes  and  is  divided  into  21  span-wise 
(3.0  in  apart)  stations.  Models  with  finer  meshes  were  used 
in  the  early  stages  to  check  for  convergence.  The  mesh 
selected  seemed  to  be  a  reasonable  compromise  between  the 
computational  time  and  accuracy. 

The  orthotropic  material  properties  of  the  lamina  are 
listed  below: 


En,psi . 23.1  x  106 

E22iPsi . 1.4  x  106 

V\2 . . . . . ..0.338 

G12,psi . 0.73  x  106 
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TABLE  1  -  MODEL  DEFINITION 


Example  GRIDS,  POINTS  AND  LINES 

Vector  5  L45  4  L14  1  L51  5 


+  20  k 
-70 
+  20 
-70 
-70 


G3 

G1 

G2 

L512 

L411 

L18 

L512 

12  L1112  11  L811  8  L128  12 


+  20  j 

POSSIBLE 

MODELS 


(not  to  scale,  e.g.  line  L512  is  30.0  in  length, 
L45  is  .286  in,  and  L14  is  about  1.15  in) 

*  =  vector  inward  ;  o  =  vector  outward. 

RIGHT  HAND  RULE  NORMAL  VECTOR  TO  SURFACE 
GRID  DEFINITION 


Ml$$ 

VECTOR 

* 

* 

* 

GRID 

DEF. 

L512 , L1112 , 
L411.L45// 

L411 , L811 , L18 , 
L14/ / 

L18 , L128 , L512 , 
L51// 

M2$$ 

VECTOR 

* 

.  .  .  . . _  .  .  .. 

* 

o 

GRID 

DEF. 

L512,L1112, 

L411.L45// 

L411 ,L811 ,L18 , 
L14// 

L512 , L128 , L18 , 
L51// 

M3$$ 

VECTOR 

* 

o 

o 

GRID 

DEF. 

L512 , L1112 , 
L411.L45// 

L18 , L81 1 , L411 , 
L14// 

L512 , L128 , L18 , 
L51// 

M4$$ 

VECTOR 

o 

o 

o 

GRID 

DEF. 

L411 , L1112 , 
L512.L45// 

L18 ,L811 , L411 , 
L14// 

L512 , L128 , L18 , 
L51// 

M5$$ 

VECTOR 

o 

o  i 

* 

GRID 

L4 1 1 , Lll 1 2  , 
L51 2 , L45// 

L18 , L811 , L41 1 , 

L18 , LI  28 , L512 , 

DEF. 

LI  Ml 

L51// 

M6$$ 

VECTOR 

o 

* 

* 

GRID 

DEF. 


L41 1 , LI 1 1 2 , 
L512.L45// 


L411 ,L811 ,L18, 
L14// 


L18 , L128 , L5 12 , 
L51  // 


IV 


STATIC  AND  DYNAMIC 


ANALYS I S 

The  six  models  described  in  chapter  III  are  basically 
subjected  to  three  main  types  of  independent  static  loads. 
Further,  a  modal  analysis  is  also  performed  to  obtain  the 
first  four  fundamental  modes  of  vibration  frequencies  and  mode 
shapes . 

The  first  load  applied  is  a  torque  at  the  tip,  not 
combined  with  any  other  load. 

The  second  load  case  consists  of  combined  lift  and  blade 
weight  at  zero  RPM,  i.e.  with  no  centrifugal  forces. 

The  third  load  case  is  the  centrifugal  force  due  to 
rotation  of  the  blade.  The  lift  and  blade  weight  are  present 
and  the  RPM  varies  from  130  rad/s  to  300  rad/s. 

A.  TORSIONAL  LOAD 

A  moment  of  26.4  in-lb  is  applied  distributed  among  the 
nodes  at  the  tip  as  concentrated  moments. 

This  load  case  reveals  the  influence  of  the  laminate 
principal  orientation  angle  (LPO).  and  the  torsional  stiffness 
for  each  model.  The  different  response  of  the  six  different 
models  also  is  presented  for  comparison  purposes. 
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The  results  are  presented  in  Figure  15  for  the  six  models 
for  LPO’s  of  +20°  and  -20°  in  each  model. 

As  much  as  40%  variation  in  the  response  may  be  observed 
in  the  comparisons. 


Figure  15  Angular  Displacement  at  the  Tip  Due  to  Pure 
Torque. 


B  .  LIFT  AND  BLADE  WEIGHT  LOADS 

This  loading  is  designed  to  get  the  combined  effect  of  the 
lift  and  blade  weight  given  by  Hodges f 1987 ] ,  expressed  in 
terms  of  radial  distance  of  the  center  as  follows: 

(17)  1  =  (0.02222  lb/in)  x  r  -  0.0123  lb/in 
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This  load  is  applied  vertically  in  two  spanwise  rows  at 
the  leading  and  trailing  edges,  in  a  proportion  that  produces 
zero  moments  at  the  quarter  chord.  This  load  distribution  is 
more  homogeneous  and  induces  less  cross-section  distortion, 
as  opposed  to  the  application  along  a  line  on  upper  surface 
as  in  beam  type  elements.  This  distributed  line  loading 
produces  small  torsional  displacement  (10“4  degree). 

The  results  are  presented  in  Figure  16  for  the  six  models, 
with  2  values  of  LPO  (+20°  and  -20°)  for  each  model. 


Figure  16  Displacement  at  the  Tip  Due  to  Vertical  Loads 
Along  the  Blade. 
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C.  CENTRIFUGAL.  LOAD 

Next  loading  case  investigated  reveals  the  effect  of 
varying  centrifugal  load,  combined  with  changing  the  LPO  angle 
within  each  model.  The  response  yields  twist  angle  induced  by 
extension-twist  coupling  present  in  the  structure  due  to 
asymmetric  stacking  sequence. 

The  LPO,  referred  to  the  top  layer  angle,  varies  in  the 
range  of  +  /-  90°,  the  internal  lay-up  in  the  stacking  sequence 
remaining  constant,  as  previously  defined. 

After  some  experiments  to  localize  inflection  points  in 
the  twist  angle  variation,  the  LPO  to  be  investigated  were 
chosen  in  the  range  of  +/-  90°,  giving  a  total  of  11  laminate 
principal  orientation  angle  for  each  of  the  six  models. 


Figure  17  Overall  Range  of  the  Twist  Angle  within  the  Six 
Models. 
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The  loads  due  to  lift  and  blade  weight  are  present  but 
their  effects  are  observed  to  be  negligible  compared  with 
those  due  to  centrifugal  loads. 

The  RPM  is  varied  from  130  rad/s  to  300  rad/s  in 
increments  of  10  rad/s,  in  a  total  of  19  load  cases.  The 
combination  of  load  cases  and  LPO  give  a  total  of  201  points 
of  twist  angle  for  each  of  the  six  models.  These  points  are 
presented  in  six  3-D  graphs  to  better  appreciate  the  choices 
offered  in  this  design  approach.  Those  graphs  present  the 
trends  of  the  twist  angle  with  RPM  and  LPO  variations.  The 
twist  angle  range  for  each  model  is  in  Figure  17,  while  the 
overall  parametric  information  are  presented  in  Figures  18 


through  23. 


Figure  18  Model  One,  Twist  Angle  Induced  by  Rotational 
Inertia  Loads  Combined  With  Laminate  Principal  Orientation 
( LPO ) ,  Referred  to  Top  Layer  Angle. 
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Figure  19  Model  Two,  Twist  Angle  Induced  by  Rotational 
Inertia  Loads  Combined  With  Laminate  Principal  Orientation 
(LPO),  Referred  to  Top  Layer  Angle . 


Figure  20  Model  Three,  Twist  Angle  Induced  by  Rotational 
Inertia  Loads  Combined  With  Laminate  Principal  Orientation 
(LPO),  Referred  to  Top  Layer  Angle, 
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Figure  21  Model  Four,  Twist  Angle  Induced  by  Rotational 
Inertia  Loads  Combined  With  Laminate  Principal  Orientation 
(LPO),  Referred  to  Top  Layer  Angle. 


k 
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Figure  22  Model  Five,  Twist  Angle  Induced  by  Rotational 
Inertia  Loads  Combined  With  Laminate  Principal  Orientation 
(LPO) ,  Referred  to  Top  Layer  Angle. 
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Figure  23  Model  Six,  Twist  Angle  Induced  by  Rotational 
Inertia  Loads  Combined  With  Laminate  Principal  Orientation 
(LPO),  Referred  to  Top  Layer  Angle. 
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D 


MODAL  ANALYSIS 


To  observe  the  effects  of  each  configuration  (Model)  and 
laminate  principal  orientation  on  the  first  four  fundamental 
frequencies,  a  free  vibration  analysis  is  performed. 

The  results  of  the  modal8  analysis  are  depicted  in  Figures 
24  and  25  for  the  six  models  and  for  LPO’s  of  -20°  and  +20° 
respectively . 


Figure  24  Frequency  Variation  among  Different  Models. 


6  1  FW  means  First  Flapwise;  1  CW  means  First  Chordwise 
and  so  on. 


50 


It  may  be  noted  that,  even  for  simple  models,  this 
solution  procedure  is  a  time  consuming  task,  only  recommended 
for  final  or  quasi  final  configurations  in  the  design  process, 
or  using  a  very  crude  mesh  lay-out. 

It  can  be  seen  that,  for  each  vibration  mode,  the 
frequency  remains  practically  constant.  This  reveals  a  low 
sensitivity  of  the  frequency  to  changes  in  model 
configurations,  considering  the  first  four  fundamental  modes 
and  same  LPO’s. 
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V 


ANALYSIS  OF  RESULTS 


In  this  chapter  the  results  obtained  for  each  model  will 
be  presented  and  several  observations  are  inferred. 

The  present  research  was  conducted  to  get  an  insight  of 
the  trends  implicated  into  the  two  questions  posed  in  Chapter 
III,  sections  B  and  C.  This  study  also  provides  some  insight 
into  answering  the  questions  posed  in  sections  B  and  C  of 
Chapter  III. 


What  is  the  magnitude  of  twist  available  from  a  particular 
geometry  and  material  system? 


What  are  the  possibilities  to  model  different  structural 
configurations  starting  with  the  same  number  of  layers, 
same  lay-up  design,  same  thickness  and  material 
properties? 


The  first  question  is  a  design  approach  to  structural 
analysis.  That  approach  will  be  elaborated  in  this  chapter. 
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The  second  question  opened  a  vast  field  for  research 
activity  that  will  be  addressed  in  next  Chapter. 

It  can  be  seen  in  Figure  14  that  the  models  1  through  6 
are  obtained  simply  by  changing  one  of  the  normal  vectors 
associated  with  the  laminate,  (in  one  of  the  three  faces  of 
the  spar).  The  effects  of  these  changes  in  the  structural 
behavior  of  the  blade  may  be  divided  into  two  main  categories, 
static  and  dynamic  behavior. 

A.  DYNAMIC  BEHAVIOR 

The  free  vibration  characteristics  are  shown  in  Figures 
24  and  25.  Figure  24  shows  the  frequencies  for  LPO  of  -20°  and 
Figure  25  presents  for  LPO  of  20°.  The  first  four  fundamental 
modes  of  vibration  have  very  little  sensitivity  to  the  model 
type,  and  may  be  observed  that  the  frequencies  do  not 
change,  for  the  same  or  symmetric  LPO  (positive/negative). 

This  behavior  can  be  considered  benign,  or  in  other  words 
auspicious,  in  a  connotation  that,  once  the  designer  arrived 
at  the  desired  frequency  values  for  the  RPM  range  considered, 
the  LPO  or  the  model  can  be  changed  until  a  desired  twist- 
extension  coupling  is  achieved. 

These  iterations  can  be  done  without  looking  very  close 
to  resonance  problems,  which  must  be  checked  anyway  when  the 
final  model  is  to  be  selected. 
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This  assumption  is  made  based  on  the  validity  of 
Yntema[1955]  approach  for  these  models.  It  is  possible  that 
the  change  in  the  stiffness  matrix  due  to  RPM  affects  the 
rotor  built  of  laminated  composite  materials  in  a  different 
way  than  isotropic  rotors.  Further  research  needs  to  be  done 
along  this  direction  to  obtain  comparable  results  for 
composite  rotors. 

In  Figures  24  or  25,  each  column  corresponds  to  one  mode 
of  vibration  at  zero  RPM  and  is  a  point  in  the  3-D  plot  for 
that  specific  mode  (flatwise,  chordwise  or  pitchwise9).  The 
RPM  and  LPO  may  be  varied  as  in  the  static  case  yielding 
carpet  plots.  The  resulting  plots  would  validate  or  not  the 
assumption  of  Yntema’s  values  for  composite  blades. 

B.  STATIC  BEHAVIOR 


The  static  behavior  is  presented  in  Figures  18  through  23. 
The  parameters  considered  are  RPM  and  LPO  as  inputs  and  the 
twist-extension  coupling  manifested  through  the  twist  angle 
at  the  tip,  as  the  output.  This  angle  will  be  referred  as  TTA 
(tip  twist  angle)  in  the  subsequent  discussion. 

The  TTA  can  be  used  as  a  measure  of  the  twist  angle  ratio 
along  the  blade  span.  The  twist  angle  is  not  a  linear  function 
of  the  radius  but  may  be  represented  by  a  second  degree 

9  Pitchwise  mode  is  a  higher  mode  in  this  case  and  is  not 
considered  in  this  report. 
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polynomial,  such  that  the  TTA  is  a  good  parametric  factor  to 
infer  the  twist  angle  magnitude  along  the  blade. 

As  mentioned  before,  the  six  models  are  obtained  by 
changing  the  orientation  of  the  normal  vector  associated  to 
the  laminate,  in  the  three  shells  that  constitute  the  spar, 
one  at  a  time  as  shown  in  Figure  14. 

In  that  figure  it  can  be  verified  that  Model  1  and  Model 
4  can  be  built  with  only  one  unbroken  piece  of  laminate, 
wrapped  around  a  mold.  Further,  one  model  is  just  the  inverse 
of  the  other,  in  other  words  Model  1  is  the  inside  out  of 
Model  4  . 

The  models  2  and  5  can  also  be  considered  as  some  sort  of 
inverse,  in  a  broader  sense;  similar  reasoning  may  be  extended 
to  models  3  and  6. 

Kith  the  above  considerations,  the  six  models  can  be 
analyzed  in  pairs,  which  is  very  convenient  from  the  point  of 
view  of  getting  the  information  in  a  simple  and  concise  form, 
so  that  the  designer  can  quickly  evaluate  the  model,  keeping 
in  view  the  design  requirements. 

The  second  important  consequence  of,  this  inverse  concept 
is  that ,  Model  1  and  Model  4  can  supply  a  spectrum  of  twist 
angle  requirements,  for  a  given  loading  condition. 

The  importance  of  this  reside  in  the  manufacturing 
process,  in  that  both  models  may  be  built  by  filament  wound 
technique.  This  process  is  less  expensive  and  more  reliable 


as  it  eliminates  the  filament  discontinuity  in  the  shell’s 
junction,  and  requires  less  accurate  quality  control. 

The  convention  used  in  the  modeling  process  has  as  a 
aftermath  that  positive  angles  have  active  aerodynamic 
stabilizing  effect,  acting  as  a  propeller  governor  device. 

1 .  Models  1  and  4 

Figures  18  and  24  exhibit  the  behavior  of  these  two 

models . 

The  models,  which  were  referred  as  inverse  to  each  other, 
also  show  the  response  characteristics  and  behave  as  expected. 

This  may  be  noted  by  observing  the  displacement  for  a 
given  RPM  and  Layer  Principal  Orientation  Angle  in  Figure  26. 

2.  Models  2  and  5 

Those  models  have  the  same  inverse  relationship  in 
the  TTA  value.  The  main  difference  is  that  these  two  models 
are  stiffer  than  the  other  four,  as  shown  in  Figure  17.  Figure 
27  shows  these  two  models  drawn  together. 

3.  Models  3  and  6 

The  inverse  relationship  also  holds  here  and  the 
stiffness  values  are  higher  than  that  found  in  models  i  and 
4  and  lower  than  models  2  and  5  (Figure  17).  See  Figure  28 
for  this  comparison. 
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Figure  26  Models  1  and  4  static  behavior  comparison ,  inputs 
are  RPM  and  Laminate  Principal  Orientation  (LPO)  angle, 
output  is  Tip  Twist  Angle  ( TTA ) . 
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Figure  27  Models  2  and  5  static  behavior  comparison ,  inputs 
are  RPM  and  Laminate  Principal  Orientation  ( LPO )  angle, 
output  is  Tip  Twist  Angle  (TTA). 
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Figure  28  Models  3  and  6  static  behavior  comparison ,  inputs 
are  RPM  and  Laminate  Principal  Orientation  (LPO)  angle, 
output  is  Tip  Twist  Angle  (TTA). 
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VI 


CONCLUSIONS  AND  SCOPE 


FOR  FUTURE  RESEARCH 

A.  CONCLUSIONS 

This  report  presents  a  detailed  analysis  of  a  typical 
rotor  blade,  the  focus  being  the  static  behavior  in  the 
presence  of  extension-twist  coupling  due  to  asymmetric 
stacking  sequence  of  laminas. 

Different  extension-twist  coupling  effects  are  obtained, 
by  changing  the  ply  orientation  with  respect  to  the  body  axes, 
while  keeping  the  stacking  sequence,  thickness  and  internal 
lay-up  sequence  as  invariants. 

These  effects  are  presented  in  terms  of  carpet  plots.  The 
extension-twist  coupling  is  measured  through  the  angular 
displacement  at  the  tip  of  the  blade  as  function  of  RPM  and 
Laminate  Principal  Orientation  (  LPO  )  angle  variations. 

Free  vibration  analysis  is  performed  for  two  LPO’s  of  the 
six  different  models.  The  first  four  fundamental  frequencies 
show  small  sensitivity  to  these  specific  structural  changes. 

The  static  analysis  revealed  a  sort  of  "antisymmetric" 
behavior  within  the  models,  allowing  them  to  be  grouped  in 
pairs.  Such  behavior  can  be  of  use  in  the  design  process,  to 
trim  structural  response  to  prescribed  loads. 
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A  structural  analyst  should  be  very  careful  in  dealing 
with  composite  laminated  structures,  especially  if  the  analyst 
is  experienced  only  in  structural  analysis  of  isotropic 
materials.  It  was  noted  that  it  is  very  easy  to  be  mistaken, 
mainly  if  there  are  asymmetries  present  within  the  laminate. 
It  is  very  important  to  know  how  the  local  coordinates  systems 
are  established  in  the  finite  element  program  being  used. 
These  axes  will  determine  the  relative  position  of  the 
laminate  with  respect  to  the  structure.  If  asymmetries  are 
present,  it  is  possible  to  get  different  structural  response 
and  have  different  structural  configurations,  all  starting 
with  the  same  basic  laminate. 

B .  SCOPE  FOR  FUTURE  RESEARCH 

Future  research  may  be  beneficial  into  two  main  areas. 

The  first  one  being  to  validate  these  models  with 
experimental  data,  providing  the  necessary  confidence  in 
further  utilization  of  the  method  envisaged. 

Once  the  validation  is  done,  the  method  can  be  corrected 
and  trimmed  for  operational  use. 

Then  a  field  that  calls  for  attention  is  a  dynamic 
characteristics  mapping,  that  is,  with  the  RPM  and  LPO  as 
parameters  obtain  the  fundamental  frequencies  as  output. 

This  would  be  a  3-D  composite  material  extension  ot 
Yntema’s  [1955]  work  for  beams  of  isotropic  material. 
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ANNEX  1 


Ml  2  O  .  SRC 


£  DEFINES  BLADE  MODEL  WITH  NO  ELIPSE  PROJECTION 

£  (  LOST  3.976  DIGITS  IN  DECOM  ) 

£  (  ~  20  MIN  TO  RUN  2  LDCASES  ) 

£  NO  ROOT  DEFINED  JUST  SUPRESSED  DOF  OF  FIRST  RIB 

£ 

KPOINT 
1/. 65, 5. 23,0 
2/ .455,5.23, . 1103 
3/0,5.23, . 1545 
4/-. 49, 5. 23, .1433 
5/-. 49, 5. 23,-. 1433 
6/0,5.23,-. 1545 
7/. 455, 5. 23, -.11 03 
8/. 65, 35. 23,0 
9/. 455, 35. 23, . 1 103 
10/0,35.23, . 1545 
11/-. 49, 35. 23, . 1433 
12/-. 49, 35. 23,-. 1433 
13/0,35.23,-. 1545 
14/. 455, 35. 23, -.1103// 

$ 

£  LINES  GENERATION 

£ 

SLINE 
L34/3 , 4 , 2 
L45/4 ,5,4 
L56/5 , 6 , 2 
LI  8/1 ,8,21 
L41 1/4, 11,21 
L512/5, 12,21 
L101 1/10, 1 1 ,2 
LI  1 12/1 1 , 12,4 
L 1 21 3/ 1 2 , 13,2// 

£ 

£  ARC  GENERATION 

£ 

CARC 

C13/1 ,2,3,5 
C6 1/6,7, 1  ,5 
C81 0/8 , 9 , 10,5 
0138/13,14,8,5// 

£ 

£  COMPOSITE  LINES  GENERATION 

£ 

COMPLINE 
L 14/C  1 3 , L34 
L51/L56,C61 
LSI i/C8!0,L!01 1 
L 1 26/ L 121 3 , C138// 


*  BLADE  MATERIAL  AND  ELEMENTS  DEFINITION  (IM6/R6376) 

$  BLADE  MASS  DISTRIBUTION  =  .0123  LB/IN  =  .00038233  SLUG/IN 

$  (COMPOSITE  DENSITY: .001796  SLUG/INA3) 

S  (TUNGSTEN  DENSITY  =.021116  SLUG/INA3) 

*  (MEAN  DENSITY  =.004418  SLUG/INA3) 

*  (COMPOSITE  AREA  =.086525  INA2) 

$  (COMPOSITE  VOLUME  =2.5957  INA3) 

$  (TUNGSTEN  AREA  =.01358  INA2  =>  DIA=.1315  IN) 

*  (MODEL  FOCUSED  LENGHT  =  30  IN;  TOTAL  LENGHT  =  35.23  IN) 

ORMAT , 53 

1/5 . 07E5 , 23 . 1 E6 , .338, .73E6, .004418 
8. 12E3, 1 .4E6, 1 .41E4// 

$ 

$  LAYER  DEFINITION  O20/-70/+20/-70/-70/+20] 

$  THICKNESS  =  .0055  IN/LAYER  *  6  =  .033  IN 

COMPTH 

1,6/1, 1,1, 1,1,1 
.0055, .01375,20 
.0055, .00825,-70 
.0055, .00275,20 
.0055, -.00275, -70 
.0055, - .00825, -70 
.0055, -.01375,20// 

S 

$  ELEMENTS 

*  GRIDS  MUST  BE  GENERATED  IN  SAME  ORIENTATION 

if  STABLISW  A  VIEW  POINT  NORMAL  TO  SURFACE  OUTSIDE  (OR  INSIDE) 

*  WITHOUT  CHANGING  THIS  REFERENCE  GENERATE  THE  GRIDS  FOR 

$  ORTHOTROPIC  MATERIALS  (FILAMENT  WOUNDED  THINKING) 

S 

$  THE  FIRST  SIDE  OF  GRID  DETERMINES  THE  REFERENCE  FOR  MATANG 

$ 

GETY/TB3/1 , 1// 

GR I D4 

G1/L41 1 ,LS1 1 , L IB, L 1 4 
G2/ L 1 8 , L 1 28 , L5 1 2 , L5 1 
G3/L5 1 2,L1112,L411 ,L45// 

END 


ANNEX  2 


LAYUE . AUX 


FILE  TO  BE  USED  IN  BUILDING  THE  SEVERAL  LAYUPS  IN 
EACH  M$**.SRC  FILES 
$  =  LAYUP  NUMBER 

**  =  ANGLE  OF  TOP  LAYER  WITH  RESPECT  TO  NORMAL  VETOR  DIRECTION 

1-  LAYUP  GRID  DEFINITION  --  6  MODELS  -- 

Ml**. SRC 
GRID4 

G1/L411,L811,L18,L14 
G2/L 18,L128,L512,L51 
G3/L51 2 ,L1112,L411 ,L45// 

M2**. SRC 
GRID4 

G1/L41 1 ,L81 1 , L 18, L 14 
G2/L51 2 ,L128,L18,L51 
G3/L51 2 ,L1112,L411 ,L45// 

M3**. SRC 
GRID4 

G1/L18,L811,L411,L14 
G2/L51 2 ,L128,L18,L51 
G3/L5 1 2 ,  L  1 1 1 2  ,  L4 1 1 ,  L45//' 

M4** . SRC 
GRID4 

G1/L1S,LS1 1 , L  4 1 1 ,L14 
G2/L512 , L128,L18,L51 
G3/L411 , L 1 1 1 2 , L51 2 , L45// 

M5** . SRC 
GRID4 

G1/L18,L8!1 , L  4 1 1 ,L14 
G2/L 18,L128,L512,L51 
G3/L41 1 , L 1 1 1 2 , L51 2 , L45// 

M6** . SRC 
GRID4 

G1/L411 , L8  3 1 ,L18,L14 
62/  L  1 S  ,  L  1 26 ,  L  5  3. 2  ,  L  5 1 
53/ L 4  I  1  ,  L  1  1  12 ,  LlTl2 ,  L45// 


2-  TOP  LAYER  ANGLE  --  18  ANGLES  -- 

******************************************** 
M$00 . SRC 
.0055, .01375,0 
.0055, .00825,-90 
.0055, .00275,0 
.0055, -.00275, -90 
.0055, - .00825,-90 
. 0055 , - .01375 , 0// 

***************-r**************************** 

MS10.SRC 
.0055, .01375, 10 
.0055, .00825,-80 
.0055, .00275, 10 
.0055, -.00275, -80 
.0055, -.00825, -80 
.0055, - .01375, 10// 

******************************************** 

M$20 . SRC 

.0055, .01375,20 

.0055, .00825, -70 

.0055, .00275,20 

.0055,  -.002*75,  -70 

.0055, -.00825, -70 

.0055, -.01375,20/7 

*  ******************************************* 

M* 30. SRC 

.0055, .01375,30 

. 0055 , . 00825 . -60 

.0055, .00275,30 

. 0055 , - . 00275 , -60 

. 0055 , - . 00825 , -60 

. 0055 ,-.01 375 , 30/ / 

******************************************** 

M$40 . SRC 

.0055, .01375,40 

.0055, .00825,-50 

. 0055 , . 00275 , 40 

. 0055 , - . 00275 , -50 

. 0055 , - . 00825 , - 50 

. 0055 ,-.01 375 , 40// 


M$50„SRC 
.0055, .01375,50 
.0055, .00825,-40 
.0055, .00275,50 
.0055, -.00275, -40 
.0055, -.00825, -40 
.0055, -.01375,50// 

##*###  *************************************** 

M$60. SRC 

.0055, .01375,60 

.0055, .00825,-30 

.0055, .00275,60 

.0055, -.00275, -30 

.0055, -.00825, -30 

.0055, -.01375,60// 

#***#***#******##*#*#•#••*•****#***#****•****.*#** 

MS70.SRC 
.0055, .01375,70 
.0055, .00825,-20 
.0055, .00275,70 
.0055, -.00275,-20 
.0055, -.00825, -20 
.0055, - .01375,70// 

******************************************** 

M$80 . SRC 
.  OC'SS  ..01  375 , 80 
.00825,-10 
.0055,.  "75,80 

.0055,-  J275 , - 1 0 

.0055, -.00825, -10 
.0055, -.01375,80// 

******■»**•»****#*!>  •  *■#********.»****♦•****_#*_***•» 

Mi90.SRC 
.0055, .01375,90 
.0055, .00825,0 
.0055, .00275,90 
. 0055 , - . 00275 , 0 
.0055, -.00825,0 
•  0055,-.013*’5,90// 


**********#*+**#********♦****#*_**.#********:** 

Mi -90. SRC 
.0055, .01375,-90 
.0055, .00825,0 
.0055, .00275,-9; 

.0055. - . 20275,0 
. 0055 , - . 00825,0 
.  ‘*055  ,-.01 375 ,  -*50/  . 


7  U 


MS-BO. SRC 
.0055, .01375,-80 
.0055, .00825.10 
.0055, .00275,-80 
.0055, -.00275, 10 
.0055, -.00825, 10 
.0055,- .01375,-80// 

******************************************** 

MS-70. SRC 
.0055, .01375,-70 
.0055, .00825,20 
.0055, .00275,-70 
.0055, -.00275, 20 
. 0055 , - . 00e25 , 20 
.0055, -.01375, -70// 

************************* ******************* 

MS-60 . SRC 
.0055, .01375,-60 
.0055, .00825,30 
.0055, .00275,-60 
.0055, -.00275, 30 
.0055. -.00825, 30 
.0055, -.01 375, -60// 

*********************#************-********** 

MS’ -50.  SRC 
.0055, .01375,-50 
.0055, .00825,40 
.0055, .00275,-50 
.0055,- .00275,40 
.0055, -.00825, 40 
.0055, -.01375, -50// 

******************************************** 

MS-40. SRC 
.0055, .01375,-40 
.0055,-00825,50 
.0055, .00275,-40 
.0055, - .00275,50 
.0055, -.00825, 50 
.0055, -.01375, -40// 

*******#***********************:******:******# 

MS-30, SRC 
.0055 ,.01375,-30 
.0055. .00835,60 
.0055, .00275,-30 
.0055, - .00275,60 
.  0055 ,  -  .  00825 ,  t>0 


******************************************** 

MS-20 . SRC 

.0055, .01375,-20 

.0055, .00825,70 

.0055, .00275,-20 

.0055, - .00275,70 

.0055, - .00825,70 

.0055, -.01375, -20// 

******************************************** 
MS- 10 . SRC 

.0055, .01375,-10 
.0055, .00825,80 
.0055, .00275,-10 
.0055, - .00275,80 
.0055, -.00825, 80 
.0055, -.01375, -10// 
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